Topological regulation of activation barriers on fractal substrates 
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We study phase-ordering dynamics of a ferromagnetic system with a scalar order-parameter on 
fractal graphs. We propose a scaling approach, inspired by renormalization group ideas, where a 
\ crossover between distinct dynamical behaviors is induced by the presence of a length A associated to 

, the topological properties of the graph. The transition between the early and the asymptotic stage 

D ■ is observed when the typical size L(t) of the growing ordered domains reaches the crossover length 

A. We consider two classes of inhomogeneous substrates, with different activated processes, where 
the effects of the free energy barriers can be analytically controlled during the evolution. On finitely 
ramified graphs the free energy barriers encountered by domains walls grow logarithmically with L(t) 
while they increase as a power-law on all the other structures. This produces different asymptotic 
growth laws (power-laws vs logarithmic) and different dependence of the crossover length A on the 
model parameters. Our theoretical picture agrees very well with extensive numerical simulations. 
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I. INTRODUCTION 
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The slow relaxation of glasses, spin-glasses and phase-separating systems is a subject of paramount importance in 
, non-equilibrium statistical mechanics. In the case of homogeneous ferromagnets, i.e non-disordered magnetic systems 

defined on homogenous substrates such as lattices, the basic features of the dynamics after a quench below the 
critical temperature are well understood in terms of the domain growth mechanism The unbounded growth of 
^ • the size of ordered domains with nearly equilibrium composition involve a scaling symmetry due to the presence of 
a dominant length scale, which is manifested in the lack of time-translational invariance and aging. Although this 
simple paradigm of slow relaxation is expected to encompass a broad variety of situations, its applicability to more 
complex systems, such as spin-glasses, where both homogeneity and ferromagnetism are spoiled, remains a debated 
issue • An intermediate class of systems is that of non- homogeneous ferromagnetic materials Q . In these systems 
space homogeneity is wrecked by spacial modulations of some relevant parameter which, in order to maintain the 
system ferromagnetic, must coexist with the low-temperature order. In this case phase-ordering kinetics is preserved, 
although interesting new features may arise. The agent whereby homogeneity is spoiled may be random, like random 
fields, random bonds, or dilution, or deterministic, as in the case of models defined on deterministic graphs. The 
kinetics of these systems have been an active area of research for quite some time now. An unifying feature is the 
ubiquitous appearance of energetic barriers slowing down the motion of domains boundaries. This happens because, 
5J] , due to space inhomogeneity, interfaces are pinned in particular positions, and this may have dramatic consequences 
on the growth law and on the properties of correlation functions 043- Despite many experimental and theoretical 
studies |3[, a number of issues are still open. In particular, the conditions under which inhomogeneities may radically 
modify the asymptotic dynamics are not a priori known [8|. Indeed there are cases where disorder merely change 
the time units, like the one-dimensional Ising model with random bond @, and others, like the same model but with 
disorder in the form of a random-field [lOj], where profound qualitative modifications occur. For this second class of 
systems the nature of the dynamical scaling symmetry and of the asymptotic growth law is not yet fully characterized. 
Related to that, the conjecture of a superuniversal behavior, namely the idea that scaling functions are robust with 



respect to the presence on inhomogeneities which do not to change the low temperature properties of the system [11 
has been proposed. 



Recently, in the context of disordered ferromagnets, the observed behavior has been interpreted 12[ in terms 
of a picture inspired by renormalization group ideas where disorder introduces, next to the domains size, another 
characteristic length which gives rise to a crossover pattern. In this paper we apply similar ideas to the realm of 
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ferromagnets defined on non-homogeneous physical [13J fractal graphs, which turn out to be a simple but paradigmatic 
case where energy barriers thoroughly modify the asymptotic behavior of phase ordering dynamics. Indeed, at variance 
with the case of disordered systems, the nature of the energetic barriers and their scaling properties can be reasonably 
well understood and they can be related to the topological inhomogeneity. This allows one to make precise predictions 
on the growth laws, on the scaling properties and on the crossover phenomena. These predictions conform quite well 
to the outcome of numerical simulations of the dynamics of the Ising model on graphs. The main result of this paper 
is the existence of a crossover from an early power-law behavior to a slower asymptotic growth. The latter may be 
power-law (with a smaller, temperature dependent exponent) or logarithmic according to the structure of the barriers, 
which in turn depend on the topology of the graph. In particular, we conjecture power-laws to be associated to graphs 
which sustain ferromagnetic order only at zero temperature, while logarithms are expected if an ordered phase exists. 
The scaling functions of correlation functions are sensitive to the crossover. 

This paper is organized as follows. In Sec. HU we provide an overview of domain growth laws and of the scaling 
framework proposed in [l2[ for disordered systems. In Sec. IIIII we specify and adapt the scenario to the case of 
ferromagnets on graphs, providing also predictions for several quantities such as crossover lengths and exponents. In 
Sec. IIV[ we present numerical results for the phase-ordering of the Ising model on three model fractal graphs, the 
Sicrpinski carpet (SC), the Sicrpinski gasket (SG), and the T- fractal (TF). For the first structure a low-temperature 
ordered phase exists, whereas it is not sustained in the other two. We will argue this fact to represent a basic difference 
for the dynamical evolution. The numerical results are interpreted using the scaling ideas of Sees. IIIIIIII and agree 
quite nicely with the predictions obtained there. In Sec. [V] we conclude the paper with a summary and a discussion 
on the generality of our results. 



II. GROWTH LAWS AND DYNAMICAL SCALINGS 



Domain growth kinetics is characterized by an ever increasing typical domain size L(t) after the quench of a system 
to a final temperature T in the ordered phase. When L(t) becomes the dominant length in the problem, dynamical 
scaling is observed, meaning that all other lengths can be measured in units of L(t). Such property is reflected for 
instance in the two-point/two-time order parameter correlation function Q(r;t,t w ), where r is the distance between 
the two points and t > t w , which in homogeneous systems scales as [l4[ 

G(r;t,t w ) = g[r/L(t),L(t w )/L(t)}, (1) 

Q being a scaling function. 

Generally speaking, the growth law and other properties such as Q may depend on several factors, e.g. temperature, 
conservation laws, dimensionality, order parameter symmetry, disorder and topology of the substrate. In the simplest 
cases, as in homogeneous magnets with non-conserved dynamics, the kinetics proceeds without activation events. This 
means that no free energy barriers are encountered in the evolution and, consequently, phase-ordering is observed down 
to temperature T = 0. In this case the typical domains size increases as L(t) ~ t 1 / 2 . On the other hand, the kinetics 
of other coarsening systems require thermal activation. The simplest example are perhaps homogeneous magnets with 
conserved dynamics. In this case energy barriers arise as due to the microscopic evaporation-condensation mechanism 
underlying the evolution fl5l ]. In general, barriers can have different origin and scaling properties. In non- homogeneous 
systems, like the disordered ones or those embedded on a non-homogeneous substrate, they are typically due to the 
pinning of interfaces in certain preferred positions. In a wide class of disordered systems, such as the Ising model with 
random fields [Io|, [H, EE Q3 , random bonds [lU , or random dilution Q , this pinning effect slows down asymptotically 
the growth law to a logarithmic form 

L(t)^[\n(t/r)}^, (2) 
where r and ip are model-dependent quantities. With an Arrhenius form 

t oc e £ /^ (3) 

for the time needed to escape a barrier with activation free energy £, the problem of establishing the growth law is 
directly related to the one of determining the typical height £ of barriers at a certain stage of the evolution. Indeed, 
inverting the relation one concludes that the typical activation increases with the domains size as a power-law: 

£(L) oc L*. (4) 

Eq. dH is not the only possible situation, since a different algebraic growth law 

L(t)oct 1 / c , (5) 
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with C a temperature dependent exponent, is observed asymptotically in systems denned on a class of fractal structures 
fl9l |2G ] and pre-asymptotically in random field [HJ and random bond ferromagnets [HI, [2l[ . Plugging Eq. (J5J) in Eq. 
([3]) one arrives at 

£{L) oc ML)% (6) 

for the scaling of the barriers with the domains size. A classification of phase-ordering systems according to the 
possible growth laws is made in Q. In general, however, a complete characterization of L(t) in the many varied 
instances of coarsening systems is lacking. Besides that, understanding the behavior of the scaling functions such as 
Q in Eq. (TO) is also an open problem. 

In Ref. [12J , it was proposed to unify the wide variety of behaviors observed in domain growth into a generalized 
scaling framework. The basic idea, originally conceived for disordered systems, is that the agent responsible for 
inhomogencity introduces an extra characteristic length A. In the case of the Ising model with non-conscrved dynamics 
defined on a inhomogeneous graph considered in this paper (see section IIII[) , the only parameters of the model are 
the ferromagnetic coupling constant J and T, which enter in the combination e = J/T. Then it must be A = A(e). 
The presence of A introduces a crossover phenomenon between two different dynamical behaviors when L[t) crosses 
A. For the domains size this is assumed to be described by 

L{t;e) = t 1 l*L{e/t% (7) 

where z and <p arc a growth and a crossover exponent, with a scaling function behaving as 

L{x) = const., for x — > , , 

L{x) = xW'hix- 1 /*), for noo { ' 

where x = e/P'. In systems where disorder plays the role of an irrelevant parameter <j> is positive, inhomogenities only 
affects the dynamics at early times and L(t) crosses over to the pure system behavior L(t) ~ t x / z at large times (12j ]. 
This usually happens if there is an upper limit for the height of energetic barriers. On the other hand, when barriers 
of any size can be encountered, as in the fractal models considered in this paper and others, (f> is negative and Eq. 
([8]) describes the crossover from the early stage power law 

L(t)~tV», (9) 

with a temperature independent growth exponent z, to the asymptotic form L(t) = X(e)£(t/e 1 ^) at the crossover 
length A = e 1 /^ 2 ). The functions A(e) and (.{t/t 1 ^) are in general hard to be determined, both analytically and 
numerically. However, some prediction can be made in the case of coarsening on deterministic fractal structures, as 
will be discussed in Sec. IIIII 



A. Autocorrelation function 



In order to discuss the effect of inhomogeneities on the form of the correlation functions let us consider as a 
paradigm the autocorrelation function C{t,t w ) = Q{r = 0;t, t w ). Generalizing the crossover approach Q to this 
two-time quantity, one would expect the form 



C(t,it 



= C 



L(t) A(e) 



L{t w y L{t w ) 



(10) 



where C is a scaling function. In the context of disordered ferromagnetic systems a superunivcrsality conjecture was 
proposed according to which the effect of disorder is simply accounted for by the slower growth of L(t), while scaling 
functions entering correlation functions remains unchanged with respect to the clean case. This amounts to say that, 
e.g. for the autocorrelation function, C(x,y) should not depend on the second entry. The superuniversality property 
was checked in several models, arriving at different conclusions. Indeed, while d = 1 results |9| . 1 1 O f clearly demonstrate 
the absence of superuniversality, in the cases with d > 2, there is evidence both in favor [l7l. |22| and against 0, [l2[ 
its validity, and there is presently an intense debate on the subject. 

The original formulation of superuniversality was conceived for systems where inhomogeneities arc introduced by 
disorder. The problem of determining the relevance of inhomogeneities on the scaling functions, however, may be 
posed on more general grounds and the simple systems studied in this paper may help the clarification of this issue. 
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III. COARSENING ON FRACTAL STRUCTURES 

We will consider in the following the Ising model defined by the Hamiltonian 

H[<r]=-jY,<n<ri, (11) 

where Cj = ±1 is a unitary spin and < ij > are nearest neighbors on a graph. The dynamics is introduced by randomly 
choosing a single spin and updating it with a transition rate that in the numerical simulations will be chosen in the 
Metropolis form 

w([a] —¥ [a'}) = mm{l,ex-p[-AE/(k B T)}}. (12) 

Here [a] and [a'] are the spin configurations before and after the move, and AE = H[cr'] — H\a\. Phase-ordering is 
observed after a quench from an high initial temperature (assumed to be infinite in the following) to T = or to 
a temperature below the critical one T c . In order to infer the properties and the physical content of the quantities 
introduced in the previous section, let us focus on the structures of the SG, of the TF, and of the SC [H, [24| . 
representing prototypical examples of finitely ramified (the first two) and non-finitely ramified fractals (the latter). 
Indeed a wide class of fractal structures feature topological properties similar to the SG, TF and SC [2j|. In finitely 
ramified fractals, an arbitrary large part of the structure can be disconnected by removing a finite number of cutting 
bonds. In general, they admit an ordered phase only at T = 0. In a broad sense these structures can be considered a 
non trivial topological analogous of a low dimensional homogeneous systems with d < g^, di being the lower critical 
dimension. Non-finitely ramified graphs, on the other hand, can possess a low temperature ordered phase below 
a critical temperature T c . All these fractals can be build starting from an elementary object, denoted as the first 
generation G\, obtaining then an object of generation G2 by combining G\ parts, and then proceeding recursively as 
sketched in Fig. [1] In doing that, the linear size L n of the structure at generation n grows as L n = f ■ L„_i, where / 
is constant which depends on the structure considered and, specifically, f = 2 for the SG and the TF, and / = 3 for 
the SC. 




FIG. 1: Construction of the first 3 generations of the TF (upper part), SG (central part) and of the SC (lower part). 

Since the Ising model on finitely ramified fractals is characterized by T c = 0, phase-ordering could in principle be 
truly asymptotic only after a quench to T = 0, where however the dynamics is frozen because activated moves are 
prevented. Nevertheless, as explained in [lj|, a preasymptotic coarsening stage of diverging duration is observed for 
sufficiently low temperatures, similarly to what happens for instance on the Id lattice. We will focus on this stage of 
the dynamics to study coarsening on finitely ramified structures. 

During the evolution, energetic barriers arise due to the mechanism sketched in Fig. [2] Let us consider the SG 
first, represented in the upper panel. The figure shows schematically the evolution of an interface which progressively 
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spans a part of the structure. Initially, the position of the domain wall is outside the structure represented in the 
figure, in the left corner. This means that all the spins are, say, down. As time goes on, the interface enters the graph 
by moving across the intermediate position j(™ ax *> indicated by a dotted green line (this means that spins on the left 
of the green line have been reversed up). The index n refers to the fact that the interface is currently spanning the 
n-th generation of the fractal. When the spins of the whole generation n have been reversed, the interface, depicted 
with two green arches, is located in the configuration I n mm ' on the four cutting bonds. Since the energy of a domain 
wall is J times its length, it is clear that in I n mm ' the interface has a minimum energy E n mm ^ ■ On a finitely ramified 
structure, since the number of cutting bonds is finite this quantity is bounded from above and for large n, it becomes 
?i-independent, i.e. En nm ^ = E^ mm \ For the SG, i<;( mm ) = 4 J for any n, since there are four cutting bonds for the 
interface, independently of n. Let us assume that the highest energy E^ nax ^ of the system during the above process 
was reached in the (generic) configuration j^™ 1 ^^ so that a barrier of height £„ = E n max ' — En™* has been crossed. 
Now the interface must proceed again to the right in order to reverse all the spins of the next generation n + 1, thus 
reaching the position located on the other four cutting bonds, and indicated by the couple of magenta arches (all 

the spins in the figure at this stage have then been reversed). Also this configuration has an energy equal to ij( mm ). 
The topology of the finitely ramified graph allows the system to attain this configuration by sequentially reversing 
parts of generation n of the structure. For instance, in the Sicrpinski gasket, by first reversing another triangle of 
generation n, say the lower-right one in Fig. [5] This event is analogous to the one described before. In particular, the 
interface in the intermediate position depicted with a dashed- magenta line, is analogous to the previous one 

at j^ 7111 ^ (dotted-green), except for the presence of an extra part around certain cutting bonds which, in the present 
example, are indicated with a dashed-magenta arch. For low temperatures the number of these cutting bonds tends 
to be minimized and it does not depend on the generation. Denoting with E^™^ the maximum energy reached by 

the system in the reversal of the n + 1 generation, one concludes that = E n max ^ + E c , where E c is the amount 

of extra energy due to the new part of interface (the dotted-magenta arch in the figure) . Writing E c — Jn c where n c 
is the double of the number of broken bonds associated to the energy E c , for the SG one has n c = 4. We recall that 
in general n c docs not depend on n. Then 

£ n +i — £n + Jn c , (13) 
or, equivalently, rewriting £ in terms of the size L n of the n-th generation 

£(f ■ L n ) ~ £{L n ) + Jn c . (14) 
From this relation, dropping the index n, one has 

^01ni. (15) 
Inserting this form into the Arrhenius relation ([3]) we arrive at a algebraic growth law as in Eq. © with 

(16) 

and 

k B ln / 

at low temperatures. The argument has been developed referring to the deterministic SG but it is expected to hold 
also for the TF (with n c = 2) and in general for finitely ramified deterministic and disordered structures, as for all of 
them an arbitrary large part can be disconnected by cutting a finite number of links. 

A similar argument applies to non-finitcly ramified structures, and it will be schematically illustrated for the SC 
with the help of the lower panel of Fig. [2] As for the SG, an interface enters the portion of the fractal considered from 
(say) the lower-left edge, crossing a region of generation n, passing through an high energy position depicted as an 

arch-shaped dashed green line. The evolution then proceeds towards the position I n m which in the figure is depicted 
by a continuous green line that roughly spans the diagonal of the n-th generation. Here the energy reaches a minimum 
E n mm \ because the largest hole available at generation n is crossed, thus minimizing the number of disaligned spins. 
Then the domain wall moves to the next high energy configuration, depicted as a couple of dashed arch-shaped blue 
lines, and then to the minimum energy position depicted with a continuous diagonal blue line. Eventually the 
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configurations magenta are progressively reached by the interface. Notice that these configurations have the same 
energy of the blue ones. It is clear from the figure that the energies of the blue (or magenta) configurations occupied 
at generation n + 1 are the double of those of the green ones at generation n. Hence one can write, in place of (fT4|) 



£(f-L n )=u£(L n ), (18) 

where u = 2 for the SC. Hence, in place of Eqs. (fT5)l we find 

£(L n )~bJL d n f ', (19) 

where dfi = In uj In / is in general the fractal dimension of the intersection of the fractal structure with a line, that is 
the border of the interface, corresponding to the cut set (e.g. it is dfi = In2/ln3 for the SC considered in the figure), 
and from Eq. Q we arrive at Eq. ([2]), with 

V> = dfi. (20) 

This argument, developed for the SC, is expected again to be of general validity for all non-finitely ramified structures. 

In conclusion, according to the above discussion, finitely (SG and TF) and infinitely (SC) ramified fractals represent 
the typical examples where barriers increase logarithmically and algebraically with L, Eqs. (|4l [6]) respectively, due 
to their topological properties. Correspondingly the growth law is expected to be a power law (with a temperature 
dependent exponent) or a logarithmic behavior, Eqs. (|5|2|) respectively. These behaviors should be observed when 
the energy scale k B T associated to temperature fluctuations is small with respect to the height £ of the barriers. 
Conversely, for £ <C kgT, the effect of the barriers is negligible. Interfaces diffuse freely in this case and, in analogy 
to what observed on regular lattices, one expects to observe a power law behavior as in Eq. ©, characterized by the 
temperature independent exponent which describes the displacement (x 2 ) cx t 1 ^ of a random walker on the structure. 
For a fractal graph having spectral and fractal dimension d^-df respectively, it is z = 2df/d s (26|. We remark that 
this same power-law is predicted by approximate theories [271 ] to be the correct asymptotic growth-law, whereas our 
arguments indicate that it can only be preasymptotic. 

All these behaviors can be fitted into the crossover scenario described by Eq. ([7]) provided that z = d s /(2df) and 

£(x) oc x 1 ^ , for finitely ramified graphs , > 

£{x) cx (Inx) 1 ^ , for non-finitely ramified graphs, 

with £ and ip given in Eqs. (|16I20[) . 

The crossover length A is given by the condition 

£(X) ~ k B T. (22) 

Using Eqs. (fT5j) and p^|) this implies 

{A ~ exp[fcsTln //( Jn c )] , for finitely ramified graphs 

A ~ (k B T/J) 1/dfl , for non-finitely ramified graphs 



IV. NUMERICAL SIMULATIONS 



In this section we present the results of numerical simulations of the dynamics of the Ising model after a quench 
on the SG, on the TF, and on the SC. In the simulations we will consider systems with n = 9, n = 8 and n = 6 
generations respectively. With these choices finite size effects are not observed in the range of times considered. The 
system is prepared in a completely disordered state where spins are set to er^ = ±1 randomly and independently on 
each site i, corresponding to an infinite temperature equilibrium initial state. The evolution is then implemented 
flipping randomly chosen single spins with the Metropolis rate (|12[) where T is the temperature at its quench value. 
We always set fcs = 1 and J = 1 or, stated differently, we measure temperature in units of J/ks- We consider a 
modified dynamics called no-bulk-flip, where spins which are aligned with all the nearest neighborg, namely in the 
bulk of an ordered domain, cannot flip. It was shown that this dynamics, which has been tested and used in a number 
of studies pjl [H| , does not change the large-scale properties of the system and improves the speed of the simulation 
and the quality of the results. In order to grasp at least the basic physics inspiring the no-bulk-flip rule let us recall 
that in a coarsening system all the large-scale/long-time properties are uniquely determined by the motion of the 
ordered domains walls. Well inside the domains, regions of spins correlated over a length £(T) can rapidly change sign 
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over a characteristic time £(T) Z , where z is the dynamical exponent. However, since for large times t and T < T c , £(T) 
is always negligible with respect to L(t), the flipping of these thermal islands only affect small-distances/short-timcs 
properties which are not the universal ones we are interested in. Further details on this accelerated dynamics can 
be found in pjl, (28|. For the SG, we study quenches to different final temperatures in the range T <E [0.8 — 3.4]. 
For T > 3.4 the system very rapidly reaches the disordered equilibrium state and the coarsening stage cannot be 
observed. On the other hand, for T < 0.8 the dynamics becomes exceedingly slow to obtain reliable results with our 
computational resources. Similarly, for the TF we consider quenches in range of final temperatures T £ [0.2 — 3] and 
for the SCTe [1-3]. 

On a homogeneous structure, assuming scaling ([1]), the equal time correlation function 

G(r;t) = (a i (t)a j (t)), (24) 

depends only on the distance r between i and j. From this the typical length L(t) is usually extracted e.g. as the half 
height width G(r = L(t);t) = 1/2. Although on a generic graph the notion of distance is not straightforward, one 
can reduce its definition to the usual one along certain directions, as done in [2(|. Proceeding analogously here 
we arrive at a determination of L(t). This quantity is shown in Fig. [3] for quenches to different final temperatures. 

A first trivial observation is that in the structures considered in the simulations there is a lower cutoff Lhoie on the 
size of the holes. Hence for L(t) < Lhoie the fractal nature of the graph is not revealed and we expect to observe the 
same behavior as on a homogeneous lattice. Lhoie is of the order of the size of the first generation and, in the cases 
considered here, Lhoie — 3. Very early data with L(t) < L] lo i e , therefore, do not describe the effect of the fractal 
structure and we will only refer the discussion to L(t) > Lhoie- 

For all structures considered here one sees very clearly a crossover from an early regime L(t) <C A, where curves 
for different temperatures roughly collapse, to a late stage with a strongly temperature dependent behavior. This is 
what one would expect from the picture described in the previous sections. According to the discussion of Sec. IIII1 if 
such a crossover is described by Eq. (|7|), a number of predictions can be done, that we will check in the following. 

Crossover length: 

We start from the behavior of the crossover length A given in Eq. (|2"3")l . 

For the SG, using / = 2 and n c = 4 one has ln//(Jn c ) ~ 0.17. Using Eq. (|2"3"|) one concludes that, in the range 
of temperatures explored (T G [0.8 — 3.4]), A can be varied at most by a factor 1.57. In Fig. [3] (upper left panel) we 
have indicated with a dashed horizontal line the value of A obtained from Eq. (|23[) by adjusting the proportionality 
constant in such a way that, for L(t) > A only the asymptotic power-law growth (|2"Tj) (with ( ^ z) is observed. This 
can be done rather precisely for the lowest temperatures. Since A can be varied only by the small factor 1.57 it is clear 
that, even for the highest temperature one cannot obtain a significative range Lhoie < L(t) < A to clearly determine 
the preasymptotic stage, where Eq. ([9]) should hold with a temperature- independent exponent z = 2df/al s (this law 
is represented by the dashed magenta line). Notice that at the highest temperatures L(t) grows approximatively as in 
Eq. ^ but, as will be discussed below, for the different reason that the asymptotic temperature-dependent exponent 
C approaches the value z = 2df/d s at high temperatures. A similar situation is found for the TF (upper right panel). 
Here the smaller value of n c (n c = 2) allows to vary A by slightly larger factor i.e. 2.46. On this structure, then, at 
variance with the SG one should be able to detect the preasymptotic regime Lhoie < L(t) < A at least for the larger 
temperatures. 

For the SC we have plotted in Fig. [3] (second and third panel) the crossover length A obtained from Eq. (f2"3"]l . by 
using dfi = In 2/ In 3 and adjusting the proportionality constant in such a way that the crossover from early power- 
law to asymptotic logarithmic behavior occurs around A. As shown in Fig. [3] our estimate of A fit quite nicely at a 
semi-quantitative level with the data. 

Early stage: 

Next, let us consider the first regime, for L(t) <C A, that should obey the power-law behavior @ with the tem- 
perature independent exponent z = 2df/d s . As already discussed above, this preasymptotic regime is too short to 
be studied in the SG. For the TF this should in principle be observed at least at the largest temperature T = 2.5 
However, as it can be seen in Fig. [31 the crossover appears to be very broad, preventing a clearcut evidence on the 
preasymptotic behavior also on this structure. 

On the contrary, for the SC at the highest quench temperature T = 3 this stage lasts for more than a decade. Here 
one observes the expected power-law with exponent z = 2df/d s (such law is represented with the dashed magenta 
line). This shows that, indeed, in a preasymptotic regime barriers do not play any relevant role and interfaces perform 
a random walk on the graph. Notice also that, as already observed in (l9l. |20| . an oscillatory behavior is superimposed 
on a globally increasing trend. This feature is observed in all the dynamical regimes and also in the SG and TF. 



8 



Although the very limited extent of the time that can be reached in simulations does not allow to observe more than at 
most 1-2 oscillations, thus preventing any precise analysis, a semi-quantitative inspection of the data clearly suggests 
that these might be log-time periodic. This periodicity is observed in a number of apparently different phenomena, 
ranging from fracturing of heterogeneous solids [29j , to stock market indexes (30| , from magnetic systems with lack 
of translational symmetry [3l[ , to phase separating fluids under shear [32[ . This feature is generally associated to the 
presence of a discrete scale invariance (33l | which in the present case induces a recurrent trapping of the interfaces 
when a complete generation of the fractal has been ordered. The presence of such log-periodic oscilations has been 
analytically proven for a random walker diffusing on these fractal structures 34 1 . 



Late stage: 

The next step is the determination of the asymptotic growth law, namely the function £, which should behave as 
in Eq. (|2"Tj) . 

For the SG and the TF it is clear that the curves oscillate around a net power-law behavior, as already observed 
in This is what was expected for finitely ramified graphs, according to Eq. (|2"Tj) . Our arguments provide 

the prediction (|16|17j) for the low-temperature behavior of £. In order to check this we have plotted \/Q against 
temperature in the inset of the upper panels of Fig. [3] For low T the data show a good agreement with the expected 
behavior, which is represented by the dashed-blue lines. Notice that, in this case, there are no fitting parameters since 
also the value of the constant a has been inferred. 

In the case of the SC, the logarithmic growth of L(t) forces one to reach much longer times, particularly at low 
temperatures, and this in turn increases the computational effort limiting the possibility of a large statistics over 
many realizations of the thermal histories. Moreover the screening due to the oscillations is more severe. Nevertheless 
the data of the second panel of Fig [3] show quite unambiguously a crossover to a slower growth law around L(t) ~ A. 
According to our general picture this should be described by the logarithmic form of Eq. (|2~lj) . In order to check this 
we have plotted in the third panel the same data but with an extra logarithm on the time axis. In this plot the form 
of Eq. (|2"Tj) is a straight line with slope 1/ip- This is very well consistent with our data for L(t) > A. The numerical 
results show that the exponent ip decreases when T is increased. In the limit of low temperature the value of Eq. (|20|) 
is predicted. The corresponding law is represented with a dashed-magenta line in the third panel of Fig. [3] Here one 
observes a very nice agreement with the data for the lowest temperature (T = 1), suggesting the correctness of our 
argument. 



Autocorrelation function: 

Finally, our simulations allows to comment on the role of the inhomogeneities on the scaling functions of correlation 
functions. We will consider in the following the autocorrelation function 

1 N 

C{t,t w ) = —Y,{°i{t)oi{tw)), (25) 

i=l 

where N is the number of spins in the structure, for which the scaling form (|10[) is expected. 

As discussed previously, because of the smallness of A on the SG and of the broad form of the scaling functions in 
the TF, the whole crossover pattern can be better detected in the SC. Hence we concentrate on this structure in the 
following. In Fig. 0] we plot the autocorrelation function measured after different waiting times t w . The upper panel 
refers to the quench to the highest temperature T = 3, where, according to the data of Fig. |3j one can access both the 
preasymptotic regime and the crossover to the late stage, and the role of the second entry of the scaling function in Eq. 
(|10[) can be studied. Interestingly, the crossover phenomenon is fully displayed in the figure. Indeed, the two curves 
for the smallest waiting times (t w = 10 and t w = 20) almost collapse. The small residual dependence on t w , namely 
the fact that the collapse is not perfect, is probably due to the fact that the presence of the crossover is already slightly 
felt even at these early times (but we cannot reduce t w further because of the constraint L(t w ) > Lhoie)- However, the 
picture shows that the separation between these two curves is much smaller than the one among the following ones, 
and this clearly indicates a convergence of the data to a limiting mastercurve for L(t w ) <C A. In view of Eq. (fTUJ) this 
means that the second entry X/L(t w ) is so large that C [L(t)/L(t w ), X{T)/L{t w )\ ~ C [L(t)/L(t w ), oo]. Then, upon 
increasing L{t w ), the second argument of C decreases and becomes relevant. Indeed the collapse is lost and there is a 
clear tendency of the curves to move to larger values, as already noticed in fl9| . This signals the crossover from the 
early stage to the late regime. Finally, for larger values of L(t w ) (L(t w ) 3> A) one has again a tendency to collapse on 
a mastercurve which corresponds to C [L(i)/L(t w ), 0]. This is very well verified for t w = 1000, 2000, 4000. 

In order to complete the analysis we consider also the quench to the lower temperature T = 1 . In this case, since A 
is smaller the limiting curve C [L(t)/L(t w ), 0] should be achieved at earlier t w . Our results are plotted in the second 
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panel of Fig. |4j At this low temperature the autocorrelation function is strongly oscillating. This is the counterpart 
of the periodic modulations already observed in the growth-law. These oscillations hinder somehow the collapse of the 
curves. Nevertheless, one observes that already from such small waiting times as t w = 100 onwards the curves do not 
show any tendency to move upwards, at variance with the case with T = 3. This can be interpreted as due to the fact 
that, since A is much smaller at T = 1, the collapse on the mastercurve C [L(t)/ L(t w ), 0] is already achieved at these 
early times. Clearly, since A is small, the pre-crossover collapse on C [L(t)/L(t w ),oo] is not observed here. Notice 
that the two mastercurves C [L(t)/ L(t w ), oo] and C[L(t)/L(t w ),0] in the first panel are very well separated, clearly 
indicating the relevance of the second entry in the scaling function of Eq. (flQ|) . The whole behavior of C, which 
is captured by the two-parameter scaling (|10p . show unambiguously the relevance of the inhomogeneities, entering 
through the length A, in determining the shape of the scaling functions. 

V. CONCLUSIONS AND PERSPECTIVES 

In this paper we have studied the phase-ordering kinetic of a ferromagnetic system with a scalar order-parameter 
on fractal graphs. We have proposed a scaling approach, inspired to renormalization group ideas, where a crossover 
between distinct dynamical behaviors is induced by the presence of a length A introduced by the topological properties 
of the graph. The transition between the early and the asymptotic stage is observed when the typical size L(t) of 
the growing ordered domains reaches the crossover length A. In this general framework, two classes of inhomogeneous 
substrates can be defined according to the nature of the activated processes which set in during the evolution. 
Specifically, we argue that on finitely ramified graphs the free energy barriers encountered by domains walls grow 
logarithmically with L(t) while they increase as a power-law on all the other structures. This produces different 
asymptotic growth laws (power-laws vs logarithmic) and different dependence of the crossover length A on the model 
parameters. We have tested these ideas by numerical simulations of the Ising model on two model structures where, 
due to their relative simplicity, one can exhibit explicit predictions for the behavior of L(t) and of A, which conform 
very well to the numerical data. 

The models studied in this paper can be considered as simple prototypical systems to understand the more general 
and still open problem of phase-ordering in inhomogeneous systems. A natural question is then if (and how) the 
results of this Article can be extended to more general situations. For instance, one might wonder if a similar picture 
holds in systems where dilution is random instead of being, as in this paper, deterministic. Following the arguments 
developed in Sees. IIIIIIII one realizes that neither the deterministic character, neither the fractal nature are really 
determinant. Instead, the fundamental ingredient is weather the position of minimum energy of interfaces contain 
a number n c of broken bonds which is independent on L(t) or if such number scales with (some power of) L(t). 
This in turn is related to fact that the corresponding graph does not sustain a ferromagnetic phase (T c = 0) or it 
does (T c > 0), respectively. Extending this argument to the case of random dilution we can predict an asymptotic 
logarithmic growth law as in Eq. ^ for the randomly diluted (bond or site) Ising model with a fraction of occupied 
sites (or bonds) p > p c , where p c is the percolation threshold. Right at p = p c , on the other hand, we expect a 
power-law growth as in Eq. ([5])- Notice that, at variance with the fractal models studied in this Article, for randomly 
diluted systems one can tune continuously the parameter p which controls the closeness to the threshold case p = p c . 
Studies are in progress in order to test these ideas. Finally, it would be interesting to understand if similar concepts 
can be extended to interpret different system, where the inhomogeneous character is not due to dilution but to other 
agents, such as random coupling constants or spatially varying external fields. 
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FIG. 2: Schematic description of energy barriers in the SG (upper panel) and in the SC (lower panel). For the SG, the green and 
the magenta arches represent the position of the interfaces in a configuration of minimum energy after a region of generation 
n (and n + 1) has been reversed. The dashed green and magenta lines represent the configuration of the interface with larger 
energy while reversing the structure of generation n (and n + 1). For the SC, the continuous green, blue, and magenta lines 
represent the position of minimum energy of an interface after the the reversing of a part of structure of generation n (and 
?i+l). The corresponding dashed lines are the position of larger energy to be exceeded. 
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FIG. 3: The typical size L(t) (symbols) is plotted against time in a log-log plot for an Ising model quenched to different final 
temperatures (see key) on the SG (upper left panel), the TF (upper right panel) and on the SC (lower panels). The bold dashed 
magenta line in the upper and lower left panels is the short time expected behavior L(t) ~ t 1 '*, with z = 2df/d a , while in the 
lower right panel represents the asymptotic logarithmic law L(t) ~ (Ini) 1 ^, with ip = dfi = In 2/ In 3. The horizontal dashed 
lines represent the crossover length A. In the inset of the upper panels the exponent l/£, obtained from the data of the main 
part of the figure, is plotted against T. The dashed blue line is the prediction (|16I17(I . 
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FIG. 4: The autocorrelation function C(t,t m ) is plotted against L(t)/L(t w ) for different values (see key) of t w . In the first 
panel a quench to T = 3 is considered, while in the second it is T = 1. 



